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Shugo Suzuki, Shin'ichi Higai^, and Kenji Nakao 

Institute of Materials Science, University of Tsukuba, Tsukuba 305-8573 
^ National Research Institute for Metals, 1-2-1 Sengen, Tsukuba 305-0047 

(Received February 1, 2008) 

We study two-dimensional a-hole systems in boron layers by calculating the electronic struc- 
tures of Mgi_a;Naa;B2 and Mgi-^, Ala;B2. In Mgi_3;Na2,B2, it is found that the concentration 
of a holes is approximately described by (0.8 -I- 0.8a;) x 10^^ cm~"^ and the largest attainable 
concentration is about 1.6 x lO^'^ cm~^ in NaB2. In Mgi_j; Alj;B2, on the other hand, it is found 
that the concentration of a holes is approximately described by (0.8 — 1.4a;) x 10^^ cm~'^ and a 
holes are disappeared at x of about 0.6. These relations can be used for experimental studies on 
the cr-hole systems in these materials. 
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Quite recently, Nagamatsu at al. have discovered that 
magnesium diboride, MgB2, is a superconductor with a 
high transition temperature, Tc, of 39 K.^^ Extensive 
studies have now started both experimentally and the- 
oretically. In particular, since MgB2 can be regarded 
as a starting material of undiscovered high Tc supercon- 
ductors, it is important to search a variety of materials 
derived from MgB2. 

The structure of MgB2 consists of the layers of trian- 
gular lattices of Mg atoms and the layers of honeycomb 
lattices of B atoms. This is basically the same as that of 
the alkali- metal binary graphite intercalation compounds 
(GIC).^-* Since Mg and B are light elements, where the 
s and p atomic orbitals play dominant roles, the elec- 
tronic structure of MgB2^-' is also very similar to those 
of GIC.^^ In spite of these similarity, there are no GIG 
superconductors with such a high Tc] the highest Tc of 
alkali-metal binary GIG is only about 0.15 K for GgK.^^ 

One of the outstanding differences between MgB2 and 
GIG is the existence of the a holes at the center of the 
Brillouin zone,^' which are derived from the 2px and 2py 
atomic orbitals of B. Since the a bands in graphite layers 
are energetically very deep, the generation of a holes is 
extremely difficult in GIG. Furthermore, it is interesting 
to note that holes in boron layers will show character- 
istics of two-dimensional (2D) systems. As revealed so 
far, 2D systems can provide a rich variety of physics and 



possibilities of applications. It is thus important for un- 
derstanding the properties of MgB2 and its derivatives 
to study the electronic structures of the tr-hole systems 
in boron layers. 

In this Letter, we study the 2D cr-hole systems in 
Mgi_^Naj,B2 and Mgi_2,Al^B2 at x^O, 1/3, 2/3, and 
1 by calculating the electronic structures of these mate- 
rials based on the density functional theory. The main 
result of the present study is as follows. In Mgi_2;Na2:B2, 
since Na is a monovalent metal, the concentration of the 
a holes is increased with increasing x, approximately de- 
scribed by (0.8 + 0.8a;) x lO^^^m-s. in Mgi_:rAUB2, on 
the contrary, since Al is a trivalent metal, the concen- 
tration of a holes is decreased with increasing x, approx- 
imately described by (0.8 — 1.4x) x lO^^cm"'^. In the 
latter case, the a holes are disappeared at x of about 
0.6. These results can be used for experimental studies 
on the cr-hole systems in these materials. 

In the present study, we carry out first-principles cal- 
culations based on the density functional theory with 
the local density approximation, ^"^"^^ considering all elec- 
trons. To check the reliability of the results, the Kohn- 
Sham equations are solved by using both the mixed-basis 
method and the linear-combination-of-atomic-orbitals 
(LGAO) method.") In this paper, we show the results 
obtained by the mixed-basis method although the same 
results can also be obtained by the LGAO method. The 
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cut-off energy used for plane waves is 50 eV and the 

atomic orbitals employed as localized orbitals are given 
in Table I. We use not only the atomic orbitals of neu- 
tral atoms but also those of charged atoms to increase 
the variational flexibility. The number of used k points 
in the full Brillouin zone is 52 for the structure opti- 
mization of NaB2 and 185 for the electronic structure 
calculations of NaB2, MgB2, and AIB2. Also, that used 
in the calculations of Mgi_a;Naa;B2 and Mgi_2;Al2;B2 is 
104. 



Table I. Atomic orbitals used for the mixed-basis calculations. 



Atom 


Atomic orbitals (atomic charge) 


B 


ls,2s,2p(neutral) 


Na 


ls,2s,2p,3s(neutral) 


Mg 


ls,2s,2p,3s(neutral); 3p(l-|-); 3d{2+) 


Al 


ls,2s,2p,3s,3p(neutral); 3d(2-t-) 



Wc first calculate the electronic structures of NaB2, 
MgB2, and AIB2. The calculations are performed at the 
experimental lattice constants for MgB2 and AIB2; a = 
3.084 A and c = 3.522 A are used for MgB2 and a = 
3.009 A and c = 3.262 A are used for AIB2. On the other 
hand, since NaB2 is a hypothetical material at present, 
it is necessary to optimize the lattice constants of this 
material. The resultant a and c are 3.02 A and 4.19 
A, respectively, and are used for NaB2. To check the 
reliability of this result, we also optimized the structure 
of MgB2 and found that the errors for a and c are —2 % 
and -1-0.5 %, respectively. We thus believe that the result 
for NaB2 is also reliable with the same accuracy. The 
lattice constants c of these materials can be understood 
by considering the fact that the ionic radii of Na+, Mg^+, 
and Al^+ are 0.97 A, 0.65 A, and 0.50 A, respectively. 

In Figs. 1(a), 1(b), and 1(c), the calculated electronic 
structures of NaB2, MgB2; and AIB2 arc shown, respec- 
tively. The dotted lines in the figures are the Fermi level. 
The most impressive point is the position of the top of 
the (7 bands derived from the 2px and 2py atomic or- 
bitals of B. In NaB2 and MgB2, the top of the a bands 
are above the Fermi level and accordingly there exist a 
holes in these materials. This is in strong contrast to 
the fact that there are no a holes in GIG. Since Na is 
a monovalent metal while Mg is a divalent metal, the 
concentration of a holes is larger in NaB2 than in MgB2. 
Thus the concentration of the a holes can be increased 
when we increase x in Mgi_j;NaxB2. On the contrary, 
in AIB2, the top of the a bands is below the Fermi level 



(a) NaB. 
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(c) AIB2 




K r M KH A L Hr AK HM L 



Fig. 1. Band structures of (a) NaB2, (b) MgB2, and (c) AIB2. 
The dotted lines show the Fermi level. 



Two-Dimensional u-Hole Systems in Boron Layers: A First-Principles Study on Mgi_2,NaxB2 and Mgi_a:Ala!B2 



3 



and accordingly there are no u holes in this material. 
Thus the a holes are decreased when wc increase x in 
Mgi_a;Ala;B2 and eventually they are disappeared at a 
certain value of x. 

Furthermore, since the dispersion of the top of the a 
bands of all the materials is very small along the F-A 
direction, the a holes can show characteristics of 2D sys- 
tems such as large fluctuation, etc. This is in strong 
contrast to the three dimensionality of the other carri- 
ers in these materials. In all the materials, there exist 
three-dimensional (3D) tt electrons and/or holes. Also, 
in AIB2, there exist small number of 3D electrons in the 
nearly free electron state at the F point, which is derived 
from the hybridization between the 3s atomic orbitals of 
Al and the interlaycr state of boron layers; this is very 
similar to the situation in CsK, where the nearly free 
electrons also exist at the F point. -'^^^ It should be noted 
that, in GIC, there exist tt electrons and/or holes and 
also nearly free electrons and not the a holes. 




Fig. 2. In-plane (\/3 X v^) structure (dashed lines) used in the 
electronic structure calculations of Mgi_2,Ma;B2 (Af=Na or Al, 
a;=0, 1/3, 2/3, and 1) and original (1 X 1) structure (dotted 
lines). Small closed circles represent B atoms and large circles 
represent the other atoms. For a;=0, all the large circles arc 
Mg atoms. For a;=l/3, open and hatched circles are Mg and M 
atoms, respectively. For 1=2/3, open and hatched circles are M 
and Mg, respectively. For x=l, all the large circles are M atoms. 



Next, we study the electronic structures of Mgi_j,Naa;B2 
and Mgi_a;Ala;B2 at a;=0, 1/3, 2/3, and 1. In the calcu- 
lations, we assume the in-plane (v^ x v^) structure as 
shown in Fig. 2. This structure has a simplicity that the 
threefold rotation axis also exists as in the (1x1) original 



structure and thus the same Brillouin zone can be used. 

We also assumed that the lattice constants of these ma- 
terials can be obtained by linearly interpolating between 
the lattice constants of MgB2 and NaB2 or AIB2. As an 
example, we show the result for Mg2/3Na]^/3B2 in Fig. 3. 
The obtained electronic structure can be understood by 
considering the folding of the original band structures 
shown in Fig. 1. The a bands are easily identified by 
observing that they are the bands with small dispersion 
along the F-A direction just above the Fermi level. 




-16' ^ LJ I \ LJ I 1 I 1 I 1 

K r MK H A LHF AK HM L 

Fig. 3. Band structure of Mg2/3Nai/3B2. The dotted line shows 
the Fermi level. 

In Fig. 4, wc show the top of the a bands in 
Mgi_3,Naa;B2 and Mgi_a;Ala;B2 as a function of x. The 
dotted line in the figure shows the Fermi level. It is found 
that the dependence is monotonous and the largest at- 
tainable value is 1.8 eV for NaB2. We also find that, in 
Mgi_a,Ala;B2, the top of the a bands is below the Fermi 
level for x larger than about 0.6, that is, the a holes 
are disappeared for such x. The dependence on x for 
the entire region from NaB2 to AIB2 via MgB2 shown 
in Fig. 4 cannot be fit with a single straight line. It is 
necessary to fit the result with a curve or, at least, with 
two straight lines, one for Mgi_2;NaxB2 and the other for 
Mgi_xAlxB2. If we select the latter choice, the result can 
be fit with 

Etop = 0.90 + 0.91a: eV (1) 
for Mgi_a;NaxB2 and with 

£top = 0.90 - 1.57a; eV (2) 
for Mgi_3;Al3;B2. Here, we ignore the point for AIB2 in 
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obtaining the above formula because some quantities, in- 
cluding the top of the a bands and the cohesive energy 
as shown below, are not on the same straight line as the 
other Mgi_a;A^B2 are; this can be ascribed to the exis- 
tence of the nearly free electrons in AIB2 which do not 
exist in the other yigx-x^x^i calculated in the present 
study. 




_2i I I I I I I I— 

NaB2 MgB2 AIB2 



Fig. 4. Top of the u bands as a function of x. The horizontal 
axis is from NaB2 to AIB2 via MgB2. The dotted line shows the 
Fermi level. The solid line fits the four points for Mgi_j;Na3^B2 . 
The dotted-dashed line fits the three points for Mgi_3;Ala;B2, 
ignoring the point for AIB2. 



In Fig. 5, we show the cohesive energy of Mgi_a;NaxB2 
and Mgi_2;Al2;B2 as a function of x. It is found that the 
most stable material is AIB2 and the least stable one is 
NaB2. The dependence on x for the entire region from 
NaB2 to AIB2 via MgB2 can be described by a single 
straight line if we ignore the point for AIB2 because the 
cohesive energy may be affected by the existence of the 
nearly free electrons as mentioned above. The result can 
be fit with 

= 5.59 - 0.62a; eV/atom (3) 

for Mgi_3;Na2:B2 and with 

= 5.59 + 0.62a; eV/atom (4) 

for Mgi_3;Al3;B2. Although NaB2 is a hypothetical ma- 
terial at present, we believe that this material can be 
synthesized in some appropriate conditions because the 
cohesive energy of NaB2, about 5 eV/atom, is not so 
small; it is almost the same as that of the bulk Si. 
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Fig. 5. Cohesive energy as a function of x. The horizontal axis is 
from NaB2 to AIB2 via MgB2. The solid line fits the six points 
for Mgi_a,NaxB2 and Mgi_a,Ala:B2, ignoring the point for AIB2. 



Next we study the dependence of the concentration of 
a holes on a;, assuming constant density of states for the 
upper region of the a bands. This assumption is good as 
far as two conditions arc satisfied; one is that the disper- 
sion along the F-A direction is very small and the other 
is that the deviation from the in-plane free-electron-like 
dispersion is negligible. Since both conditions are sat- 
isfied as shown in Fig. 1, we derive the formula which 
give the concentration of a holes for given x. First, it is 
necessary to fit the in-plane free-electron-like dispersion 
by using effective masses for heavy and light holes. As a 
result, we find that the effective mass for heavy holes is 
0.5me and that for light holes is 0.3mc, where is the 
mass of free electrons. Next, combining with the results 
shown in Fig. 4, the following formula are derived for the 
concentration of the u holes: 

nh = (0.8 0.8a;) x 10^^ ^^-?> (5) 
for Mgi_a;Naa;B2 and 

rih = (0.8 - 1.4a;) x 10^^ ^^-i (g) 

for Mgi_a;Ala;B2. Thus, in Mgi_a;Naa;B2, the largest at- 
tainable concentration of cr holes is about 1.6 x 10^^ cm~^ 
in NaB2. In Mgi_a;Al3;B2, on the other hand, a holes are 
disappeared at x of about 0.6. 

Here, we discuss the possibility of LiB2 as a candidate 
to increase the concentration of a holes. Although one 
may expect that LiB2 is the most plausible candidate 
because of almost the same ionic radius of Li+, 0.68 A, 
as that of Mg^+ , this may not the case. We have found 
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that the structure optimization of LiB2 results in strong 

contraction of c, which is found to be less than 3 A. The 
result strongly conflicts to a simple expectation that the 
lattice constant c of LiB2 should be about 3.6 A if we 
estimate it by considering the ionic radius of Li+. To 
elucidate stable structure of LiB2, our study is now in 
progress. In spite of this result, one can still expect that 
Mgi_a;Lia;B2 for sufficiently small x can be synthesized 
because the introduction of sufficiently small Li cannot 
affect the lattice constant c so strongly. 

We next discuss the difference between MgB2 and 
other metal diborides such as transition-metal (TM) di- 
borides^^) and noble- metal diborides, AgB2 and AuB2. 
The most important point is the absence of d atomic or- 
bitals in MgB2 in contrast to the existence of d atomic 
orbitals in other metal diborides. In particular, since 
the d atomic orbitals in TM are partly filled, they form 
strong covalent bonding with a bonds of boron layers. 
This can destroy 2D cr-hole system in boron layers: wc 
have calculated the electronic structures of some TM di- 
borides and have found that the a bands of boron lay- 
ers are strongly affected by the covalent bonding with 
d atomic orbitals of TM. On the other hand, AgB2 and 
AuB2 can be candidates for similar materials as MgB2. 
The reason for this is that the a holes may survive in 
these materials because d atomic orbitals in these ma- 
terials should be sufficiently lower than the Fermi level 
and thus the hybridization between the d atomic orbitals 
and the a bands may occur at a lower-energy region. 

We finally discuss a possible relation between the su- 
perconductivity in MgB2 and the 2D cr-hole systems in 
boron layers. If the superconductivity is disappeared in 
Mgi_2;Ala;B2 at X of about 0.6, it should be caused by the 
(T-hole systems. In addition, there are no TM diborides 
with Tc as high as that of MgB2; this may be under- 
stood because they have no a holes. If this is the case, 
both of the electron-phonon interaction and the electron 
correlation can be the likely origin of the superconduc- 
tivity. Since a bonds, especially those of B, C, N, and 
O, are very strong, the interaction between the a holes 
and the in-plane cr-bond vibration is expected to be very 
strong, too. Furthermore, this electron-phonon coupling 
can result in a superconductivity of high Tc because the 
frequency of the cr-bond vibration is of the order of 0.2 
eV. On the other hand, the electron correlation in the 
(7 holes seems also important for the properties of the 
c7-hole system in MgB2 because the wave function of the 
a holes is localized to considerable degree. 
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